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We demonstrate the possibility ol a good description ol the processes <^(1020) jtvk and 
0(1020) 7/o(980) within the framework of nonrelativistic quark model assuming /o(980) to 
be dominantly quark-antiquark system. Different mechanisms of the radiative decay, that is, the 
emission of photon by the constituent quark (additive quark model) and charge-exchange current, 
are considered. We also discuss the status of the threshold theorem applied to the studied reactions, 
namely, the behavior of the decay amplitude at M^,r —> rn^ and m/g — > m^. In conclusion the 
arguments favoring the qq origin of /o(980) are hsted. 



I. INTRODUCTION 
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5^ , The if-matrix analysis of meson spectra [1-3] and meson systematics [4,5] point determinatedly to the quark- 

antiquark origin of /o(980). However there exist hypotheses where /o(980) is interpreted as four-quark state [6], 
KK molecule [7] or vacuum scalar [8]. The radiative and weak decays involving /o(980) may be decisive tool for 
understanding the nature of /o(980). 

In the present paper the reaction 0(1020) — > 7/0 (980) is considered in terms of nonrelativistic quark model assuming 
/o(980) to be dominantly the qq state, nonrelativistic quark model is a good approach for the description of the 
^ ' lowest qq states of pseudoscalar and vector nonets, so one may hope that the lowest scalar qq states are described with 
. reasonable accuracy as well. The choice of nonrelativistic approach for the analysis of the reaction (/)(1020) 7/o(980) 
was motivated by the fact that in its framework we can take account of not only the additive quark model processes 
\ (emission of the photon by constituent quark) but also those beyond it within the use of the dipole formula (the photon 
■ emission by the charge-exchange current gives such an example). The dipole formula for the radiative transition of 
[ vector state to scalar one, V — > 76* was applied before for the calculation of reactions with heavy quarks, see [9,10] 
and references therein. Still, a straightforward application of the dipole formula to the reaction 0(1020) 7/o(980) 
is hardly possible, for the /o(980) for sure cannot be represented as a stable particle: this resonance is characterized 
by two poles laying on two different sheets of the comlex-M plane , at M — 1020 — «40 MeV and M = 960 — i200 
MeV. It should be emphasized that these two poles are important for the description of /o(980). Therefore we use 
Qh! below the method as follows: we calculate the radiative transition to a stable bare /o state (this is /q '^'''^(700± 100), its 
^ ' parameters were obtained in the iiT-matrix analysis [1] . In this way we find out the description of the process 0(1020) — s- 
'T^ \ ^ fQ^^° iJOO:lsi 100) and furthermore we switch on the hadronic decays and determine the transition 0(1020) ^nn; just 
J> . the residue in the pole of this amplitude is the radiative transition amplitude 0(1020) 7/o(980). Hence we obtain 
a successful description of data for 0(1020) jmr and 0(1020) 7/o(980) within the assumption that /o(980) is 
dominated by the quark-antiquark state. 
' The conclusion about the nature of /o(980) cannot be based on the study of one reaction only but should be 
motivated by the whole aggregate of data. In the article, we list also the other processes, which provide us with 
arguments in favor of the dominant qq structure of /o(980). 

Section 2 is introductive: here we consider a simple model for the description of composite vector (V) and scalar 
{S) particles, the composite particles consisting of one-flavor quark, charge- exchange currents being absent. In such 
a model the decay transition V — > 75 is completely determined by the additive quark model process: the photon is 
emitted only by one or another constituent quark. Two alternative representations of the V — > 76* decay amplitude 
are given, namely, the standard additive quark model formula and that of a photon dipole emission, in the latter 
the factor lu ~ (my — ms) is written in the explicit form. The comparison of these two representations helps us to 
formulate the problem of application of the threshold theorem [11] to the reaction V — > 7S'. Using a simple example 
with exponetial wave functions, we demonstrate the to^ factor occurred in the partial decay width when the transition 
V — > 76" is considered in terms of the additive quark model. 

The threshold theorem has a straightforward formulation for the stable V and S states but it is not the case for 
resonances, which are the main objects of our present study. That is why we intend to reformulate the threshold 
theorem as the requirement of the amplitude analyticity - this is given in Section 3, on the basis of [12]. Working 
with nonstable particles, when V and S are resonances, the V — > 75" amplitude should be determined as a residue of 
a more general amplitude, with stable particles in the initial and final states. For example, the 0(1020) 7/0 (980) 
amplitude should be defined as a residue of the e+e~ — > 77r7r amplitude in the poles corresponding to resonances 
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^(1020) and /o(980) (the pole in the e"'"e~-channel and the /o pole in the tttt channel). The only residue in the pole 
defines the universal amplitude which does not depend on the considered reaction. 

In Section 4 we discuss the reaction </) — > 7/0 for the case, when the /o is a multicomponent system and /o and 
(j) are stable states with respect to hadronic decays. The analysis of meson spectra (e.g. see the latest if-matrix 
analyses [1,5]) definitly tells us that the /o-mesons are the mixture of the quarkonium {nn — {uu + dS)j\f2 and at) 
and gluonium components. Such a multichannel structure of /o states reveals itself in the existence of the i-channel 
charge-exchange currents. Therefore, the transition ^ 7/0 goes via two mechanisms: the photon emission by 
constituent quarks (additive quark model process) and charge-exchange current. We write down the formulae for the 
amplitudes initiated by these two mechanisms. The equality to zero of the whole amplitude at small mass difference 
of and /o, A^^^j^ ~ w at w ^ 0, is resulted from the cancellation of contributions of these two mechanisms. We 
also give a dipole representation of the transition amplitude, where A^^^^^ is determined through the mean transition 
radius and mass difference (m^ — mf^). 

Section 5 is devoted to the reaction (/>(1020) 7/o(980). First, we discuss whether it is possible to treat /o(980) as a 
stable particle. Our answer is "no" , in fact the /o(980) is unstable particle characterized by two poles located near the 
KK threshold. As was stressed above, strong transitions /o(980) — > tttt, KK reveal themselves in the two amplitude 
poles, which are located on different sheets of the complex-M plane, at M = 1020 — i 40 MeV and M = 960 — i 200 
MeV, and these two poles are important for the description of /o(980). The essential role of the second pole is seen 
by considering the tttt spectrum in ^(1020) 'ymr (Section 6): the visible width of the pick in the tttt spectrum is of 
the order of 150 MeV, and the spectrum decreases slowly with further decrease of Mt^tt- 

Another problem to be discussed is the choice of a method for the consideration of radiative decay amplitude. One 
can work within two alternative representations for the A^^^f^^ amplitude. One representation uses additive quark 
model complemented with the contribution from the charge-exchange current processes. Whence the additive quark 
model amplitude can be calculated rather definetely, at least for the lowest qq states, the charge-exchange current 
processes are vaguely determined. 

The other way to deal with the transition amplitude consists in using the dipole emission formulae, where the 
amplitude is defined by the mean transition radius and factor (m^ — mjg). For the lowest qq states, we have a good 
estimate of the radius. But there is a problem of the determination of the factor (m^ — m/p), because the /o(980) is 
unstable particle characterized by two poles. The pole M = 960 — i 200 MeV is disposed on the same sheet of the 
complex-M plane as the pole of (/)-meson, and the distance \m^ ~ ™/ol is ~200 MeV, while the pole M = 1020 — i40 
MeV is located on another sheet, the distance from the (p meson is -^70 MeV, that is also not small in the hadronic 
scale. The problem is what the mass difference factor means in case of complex masses and which pole should be 
used to characterise this mass difference. 

To succeed in the description of the decay 0(1020) 7/o(980) we use the results of the /T-matrix analysis of the 
{IJ^^ = 00^^) wave [1]. The fact is that, on the one hand, the if -matrix analysis allows us to get the experimentally 
based information on masses and full widths of resonances together with the pole residues needed for the decay 
couplings and partial widths. On the other hand, the knowledge of the if-matrix amplitude enables us to trace the 
evolution of states by switching on/off the decay channels. In such a way, one may obtain the characteristics of the 
bare states, which are predecessors of real resonances. With such characteristics, one can perform a reverse procedure: 
to retrace the transformation of the amplitude written in terms of bare states to the amplitude corresponding to the 
trasition to real resonance. Just this procedure has been applied in Section 5 for the calculation of the decay amplitude 
0(1020) -> 7/0 (980). 

Therefore, within the frame of nonrelativistic quark model, we have calculated the reaction 0(1020) 7./o"^™('^), 
where fo'^'^'^in) are bare states found in [1]. Furtheremore, with the X-matrix technique, we have taken account of 
the decays /o^'^'^(n) Tnr,KK, thus having calculated the reaction 0(1020) — > 77r7r and the amplitude of 0(1020) — > 
7/0 (980) (the pole residue in the tttt channel) . In this way we see that the main contribution is given by the transition 
0(1020) jfo'''%700 ± 100). The characteristics of /o'''''<'(700 ± 100) are fixed by the K-matrix analysis [1]: this is 
a qq state close to the flavor octet and it is just the predecessor of /o(980). In the framework of this approach we 
succeed in the description of data for the reactions 0(1020) — > 7/o(980) (Section 5) and 0(1020) jmr (Section 6). 

Let us note that such a method, the use of bare states for the calculation of meson spectra, has been applied before 
for the study weak hadronic decays £)+ — > tt+tt+tt" [13] and description of the tttt spectra in photon-photon collisions 

77 — > TTTT [14] . 

The question of what is the accuracy of the additive quark model in the description of the reactions 0(1020) 
T/o^'^C^OO ± 100) and 0(1020) ^ 7/o(980) is discussed in Section 7. We compare the results of the calculation of 
the 0(1020) 7/0 ^'^'^ (700 ± 100) reaction by using the dipole formula with that of the additive quark model. It 
is seen that, within error-bars given by the /T-matrix analysis [1], the results coincide. Still, one should emphasize 
that the dipole-calculation accuracy is low that is due to a large error in the determination of the bare-state masses. 
The coincidence of results in the dipole and additive quark model formulae should point to a small contribution of 
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processes which violate additivity such as photon emission by the charge-exchange current: this smallness is natural, 
provided the hadrons are characterized by two sizes, namely, the hadron radius {Rh ^ Rmni) and constituent-quark 
radius (r,) under the condition <^ Rf^, sec [15] and references therein. 

The performed analysis demonstrates that the studied reaction, 0(1020) 7/o(980) does not provide us any 
difficulty with the interpretation of /o(980) as Q'^-state. Still, to conclude about the content of /(980) we list in 
Section 8 the arguments in favor of the qq origin of /o(980). 



II. THE PROCESS V ^'yS WITHIN NONRELATIVISTIC ADDITIVE QUARK MODEL 

Here, in the framework of nonrelativistic quark model, we consider the transition V — > 7^ in case when the charge- 
current exchange forces are absent and the V — > -yS amplitude is given by the additive quark model contribution. 



A. Wave functions for vector and scalar composite particles 

The qq wave functions of vector {V) and scalar (S) particles are defined as follows: 



(1) 



where, by using Pauli matrices, the spin factors are singled out. The blocks dependent on the relative momentum 
squared are related to the vertices in the following way: 



2 k"^ + mev ' 
1 Gs{k^) 



2y/m fc^ -|- mev 



(2) 



Here m is the quark mass, e is the composite-system binding energy: sy = 2m — my and £5; = 2m — ms, where my 
and ms are the masses of bound states. The normalization condition for the wave functions reads 



j-^ Sp2 [*J(k)*s(k)] = J j^^Uk') SP2[(cT • fe)(eT . fe)] = 1 



d^k 
(2^ 



(3) 



Sp2 v|/+^(k)v|/^^,(k) 



d^k 
(2^ 



ipvik"^) Sp2 [a^cT^'] = S^^>. 



B. Amplitude within additive quark model 



When a photon is emitted by quark or antiquark, the V process is described by the triangle diagram, 

see Fig. la, that is actually the contribution from additive quark model. Relativistic consideration of the triangle 

diagram is presented in [16,17], while the discussion of nonrelativistic approximation is given in [12,18] (recall that in 
[18] corresponding wave hmctions were determined in another way, namely: tpvik'^) = Gv(fc^)(4fc^ + Am£v)~^ and 
^s(fc2) = Gs(fc2)(4fc2 +4mcs)"'). 

In terms of wave functions (1) the triangle-diagram contribution reads 

e^^^-^^^^r^" = eZy^,s e^Pe^S^py--^' , (4) 



^ Sp2 [*+(k)4fc„*y4k)] . 

Here elj^-* and e^a^ arc polarization vectors for V and 7: e'ii^^pv^i = and ea\a = 0- The charge factor Zv^~,'S being 
different for different reactions is specified below (see also [16,17]). The expression for the transition amplitude (4) 
can be simplified after the substitution in the integrand 
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Sp2[c7^(<7 • fe)] fc„ - ^ k^g^^, (5) 

where is the metric tensor in the space orthogonal to total momentum of the vector particle pv and photon q. 
The substition (5) results in: 

A]:^'^' = eg^^Av^,s, (6) 

where 

oo 

^Mk^)i'v{k^)^e. (7) 



The amplitudes A^^''^ and Ay^-ys were used in [16,17] for the decay amplitude 0(1020) — > 7/o(980) within rela- 
tivistic treatment of the quark transitions. 

However, for our purpose it would be suitable not to deal with Eq. (7) but use the form factor F^^''^ of Eq. (4) 
rewritten in the coordinate representation. One has: 

M/y^(k) = j d'r e*-- vl/y^(r), 

*s(k) = I d^r e*-- *s(r). (8) 

Then the form factor F^^^^ can be represented as follows: 

K^^'= l^^r Sp2 [*+(r)4fc„*y^(r)] , (9) 

where ka is the operator: fc„ = — iVa- This operator can be written as the commutator of fa and — V^/m = T 
(kinetic energy): 

2i m{Tra - r^T) = 4(-iV„). (10) 

Let us consider the case when the quark-quark interaction is rather simple, say, it depends on the relative interquark 
distance with the potential U{r). For vector and scalar composite systems we also use additional simplifying assump- 
tion: vector and scalar mesons consist of quarks of the same flavor {qq). Then we have the following Hamiltonian: 

H = - ^ + U{r), (11) 
m 

and can rewrite (10) as 

2i m{Hra - r^H) = 4(-iVa). (12) 
After substituting the commutator in (9), the transition form factor for the reaction V — > 75 reads 

F^cT''^ = J d^r Sp2 [*+(r)r„*y^(r)] 2i m{ev - eg). (13) 

Here we have used that {H + ev)^v = and (H + es)^s = 0. 

The factor ey— eg in the right-hand side (13) is a manifestation of the threshold theorem: at ey— £5) = TO5— my 
the form factor F^^^^ turns to zero. Actually, in the additive quark model the amplitude of the V — > 7>S transition, 
being determined by the process of Fig. la, cannot be zero if V and S are basic states with radial quantum number 
n = 1: in this case the wave functions ■0v'(fc^) and ipsik"^) do not change sign, and the right-hand side (7) does not 
equal zero. In order to clarify this point let us consider as an example the exponetial approximation for the wave 
functions ipvik"^) and ipsik'^)- 
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C. Basic vector and scalar qq states: the example of exponential approach to wave functions 



We parametrize the ground-state wave functions of scalar and vector particles as follows: 



1 



25/47r3/46^/* 



exp 



r 

"46^ 



(14) 



: exp 



r 



25/4773/455/4^ 

The wave functions with n = 1 have no nodes; numerical factors take account of the normalization conditions 



d\ Sp2 [M/+(r)vl/s(r)] = 1 , / Sp2 ^+^ir)^v^'{r) 



(15) 



With exponential wave functions, the matrix element for V — > -yS given by the additive quark model diagram, Eq. 
(9), is equal to 



o7/2 /,3/4l5/4 

er^4''^C""^^(additive) = (e(^)eW) 



V3 {bv + bs)'/^' 



(16) 



Formula for F^^^^ written in the frame of the dipole emission, Eq. (13), reads 

er4''^^^;cr'^"(dipole) = (e(^)e(-))^ {bf+bsfl- " 



(17) 



In case under consideration (one-flavor quarks with Hamiltonian given by Eq. (11)), the equations (16) and (17) 
coincide, F^'^T'^ (additive) = F^"^'"''^ (dipole), therefore 



m(my - ms) = by^, 



(18) 



that means that the factor es — ey in the right-hand side (13) relates to the difference between the V and S levels and 
is defined by by only. In this way, the form factor F^^^^ turns to zero only when by (or 65) tends to the infinity. 

The considered example does not mean that the threshold theorem for the reaction V ^ does not work, this 
tells us only that we should interprete and use it carefully. In the next Section, we discuss how to formulate the 
threshold theorem based on the requirement of amplitude analyticity, thus getting more information on the threshold 
theorem applicability. 



III. ANALYTICITY OF THE AMPLITUDE AND THE THRESHOLD THEOREM 

The threshold theorem can be formulated as the requirement of analyticity of the amplitude. To clarify this 
statement we consider here not only the transition of the bound states but more general process shown in Fig. 
lb, where the interacting constituents being in the vector = 1~ state emit the photon and then turn into the 
scalar = 0+ state. This amplitude has as a subprocess the bound state transition. Namely, the blocks for the 
rescattering of constituents in Fig. lb contain the poles related to bound states, see Fig. Ic, and the residues in these 
poles determine the bound-state transition amplitude (triangle diagram shown as intermediate block in Fig. Ic). 
With the notations for invariant mass squares in the initial and final states of Fig. lb as follows 

= sv, Pi = ss, (19) 

we can write the spin structures for this more general transition V — > 75. The standard representation of this 
amplitude is 

A'yr^^\sv,ss,q^ - 0) = (g^^ - ^^^) Av^-ys{sv,ss,0). (20) 

\ Sy-SsJ 
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Here we stress that the amphtude Ay^-ys describes the emission of real photon, = 0. In (20), it was taken into 
account that (Pvq) = {sv — ss)/'^- The requirement of analyticity, i.e. the absence of a pole at sy = ss, leads to the 
condition 



Av^'rs{sv,ss,0) 



0, (21) 



that is the threshold theorem for the transition amplitude V — > 'jS. 

It should be now emphasized that the form of the spin factor in Eq. (20) is not unique. Alternatively, one can 
write the spin factor as the metric tensor g^^ working in the space orthogonal to Py and q, i.e. Pv/^d'^a — ^ ^^'^ 
g^a la = 0, see Eq. (5). This metric tensor reads 

4 s 2 

5^:^(0) = 5m« + 7 '-^^—^ %<ia {Pv^qcx + QuPvoc), (22) 

[Sv - Ss)'' Sv - ss 

and we have used it in Eq. (6). The uncertainty in the choice of spin factor is due to the fact that the difference 

9^^{0) - (a^c - = 4L,40), (23) 

where 



is the nilpotent operator [12] 



S V 1 

^^"(0) = 7 <lli<la - ^Pv^tQa, (24) 

[sy-ssy 2{sv-ss) 

L^a{0)L^a{0) = 0. (25) 

The addition of the nilpotent operator _L^a (0) to spin factor of the transition amplitude V — > jS docs not change the 
expression Ay^-ysisvi ss,0), see [12] for more detail. Here, by discussing the analytical structure of the amplitude, 
it is convenient to work with the operator (20), for it is the least cumbersome. 

Consider now the reaction V ^ jS {V and S being quark-antiquark bound states), say, of the type of ^ ^ 7/0 or 
<p -yao. Because of the confinement the quarks are not the particles which form the |in) and (out| states, therefore 
the amplitudes like A^fj^^fg are to be defined as the amplitude residue for the process with the scattering of the stable 
particles, for example, for e+e~ j'K~^n~ (see Fig. 2): 

A^;:^--''-^--\sv,ss,0) = L„ - X (26) 

\ sv-ssj 

^e+e-^^ (,^_„2)(,^_^2J 5/0^.+.- +B{sy,Ss,0) 

We see that A(m|, mj^^, 0), up to the factors Gg+g-^^ and .g/„^7r+7r- 1 is the residue in the amplitude polos sy = 
and Ss = TT^'jg' just this value supplies us with the transition amplitude for the reactions with bound states cf) — > 7/0. 
If we deal with stable composite particles, in other words, if <j> and /o can be included into the set of fields |in) and 
(out I, the transition amplitude (j) 7/0 can be written in the form similar to (20); 

4t^^^''' H, ml , 0) = - ^0^7/0 H, ml , 0) , (27) 

where we have substituted Py p. For A^jt^^f^ "^/o' ^) threshold theorem is fulfilled: 

[^<P^ifo{ml,ml,0)]^, ~ {ml- ml), (28) 

4- fo 

that means that the threshold theorem of Eq. (28) reveals itself as a requirement of analyticity of the amplitude 
(j) 7/0 determined by Eq. (27). 

Let us emphasize again that formula (27) has been written for the and fo mesons assuming them to be stable, 
i.e. they can be treated as the states which belong to the sets |in) and (out|. However, by considering the process 
<p 7/0, wc deal with resonances, not stable particles, and whether this assiimption is valid for resonances is a 
question which deserves special discussion. We shall come back to this point below, and so far let us investigate how 
the requirement (28) is realized in quantum mechanics, when (f) and fo are stable particles. 
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IV. QUANTUM MECHANICS CONSIDERATION OF THE 
REACTION ^ 7F0 WITH AND Fq BEING 
STABLE PARTICLES 



In Section 2, wc have considered the model for the reaction V 76", when V and S are formed by quarks of the 
same flavor (one-channel model for V and S). The one-channel approach for ^(1020) (the dominance of ss component) 
looks acceptable, though for /o mesons it is definitely not so: scalar-isoscalar states are the mulicomponent ones. 

The existence of several components in the /o-mesons changes the situation with the 7/0 decays. First, the 
mixing of different components may result in close values of masses of the low-lying vector and scalar mesons. Second, 
equations (9) and (13) for the <j) 7/0 decay turn to be nonequivalent because of the photon emission by the i-channel 
exchange currents. 

Here we consider in detail a simple model for (f> and /o; the (p meson is treated as ss-system, with no admixture of 
the nonstrange quarkonium, nfi = {uu + dd)/^/2, nor gluonium {gg), while the /o meson is a mixture of ss and gg. 

This model can be considered as a guide for the study of the reaction 0(1020) — > 7/o(980). Indeed, the (^(1020) 
is almost pure ss state, the admixture of the nn component in 0(1020) is small, < 5%, and it can be neglected in a 
rough estimate of the 0(1020) 7/o(980) decay. 

The resonance /o(980) is a multicomponent state. Analysis of the {IJ^'^ = 00"*"^) wave in the if -matrix fit to the 
data for meson spectra tttt, KK, rjr], rjr]', tttttttt gives the following constraints for the ss, nfi and ^gf-components in 
/o(980) [1,5]: 

50% < WsslfoiQSO)] < 100%, (29) 
< Wnnif 0(980)] < 50%, 
< Wgg[/o(980)] < 25%. 

Also, the /o(980) may contain a long-range KK component, on the level of 10-20%. 

The restrictions (??) permit the variant, when the probability for the nn component is small, and /o(980) is a 
mixture of .ss and gg only. Bearing this variant in mind, we consider such two-component model for and /o 
supposing these particles are stable in respect to hadronic decays. 

It is not difficult to generalise our consideration for the three-component /o state {nn, ss and gg), corresponding 
formulae are given in this Section too. 



A. Two-component model (ss, gg) for /o and 

Now let lis discuss the model, where /o has two components only: strange quarkonium {ss in the P wave) and 
gluonium {gg in the S wave). The spin structure of the ss wave function is written in Section 2: it contains the 
factor {(T • r) in the coordinate representation. For the gg system we have Sab or, in terms of polarization vectors, the 
convolution {e[^^e'f'^). Here wc consider a simple interaction, when the potential does not depend on spin variables, 
— in this case one may forget about the vector structure of gg working as if the gluon component consists of spinless 
particles. As concern 0, it is considered as a pure ss state in the S wave, with the wave function spin factor cr^, 
see Section 2. So, the wave functions of /o and mesons are written as follows: 



The normalization condition is given by Eq. (15), with the obvious replacement: \E's ^ ^E'/p and SPy^ 
The Shrodinger equation for the two-component states, ss and gg, reads 



Us 



Ar) 



UsS^gg (r) , „ + Ugg^gg {t) 



*.5(r) 
*9p(r) 



E 



*as(r) 
*««(r) 



(31) 



Furthermore we denote the Hamiltonian in the left-hand side (31) as Hq. 

We put the gg component in to be zero. It means that the potential Uss^gg{r) satisfies the following constraints: 
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{0+ss\Us-s^gg{r)\0+gg) ^ 0, {1- ss\Us-s^gg{r)\l- gg) = 0. (32) 

These constraints do not look surprising for mesons in the region 1.0-1.5 GeV because the scalar glueball is located 
just in this mass region, while the vector one has considerably higher mass, ~2.5 GeV [19]. 

The f-exchange diagrams shown in Fig. 3a, 36,3c are the example of interaction leading to the potentials Ugs^ssir), 
Ugg^gg{r) and Uss^gg{r). The potential Uss^gg{r) contains the f-channel charge exchange. 



1 . Dipole emission of the photon 



7/0 decay 



The Hamiltonian for the interaction of electromagnetic field with two-component composite systems (quarkonium 
and gluonium components) is presented in Appendix A. 

For the transition V jS, keeping the terms proportional to the charge e, we have the following operator for the 
dipole emission: 



da 



'^'^aU ss^gg (^) ? 



, iV(xUss—^gg(j*^ 







The transition form factor is given by the formula similar to Eq. (9) for the one-channel case, it reads 



-I 



<fr Sp2 



*+(r)2rf„*^M(r) 



(33) 



(34) 



Drawing explicitly the two-component wave functions, one can rewrite Eq. (34) as follows: 



+ 



Fta^^" = /rf'r-Sp2 (r)4fc„*^(,,)^(r) 
j d^r Sp2 [*^,(gg)(r) (-ir„C/gg^,,-(r)) *^(,j)^(r) 



+ 



(35) 



The first term in the right-hand side (35), with the operator Aka, is responsible for the interaction of a photon with 
constituent quark that is the additive quark model contribution, while the term {—iraUgg^ss{r)) describes interaction 
of the photon with the charge flowing through the f-channel - this term describes the photon interaction with the 
fermion exchange current. 

Let us return to Eq. (34) and rewrite it in the form similar to (13). One can see that 



t m 



(jiofa - faHo^ = da 



(36) 



where Hq is the Hamiltonian for composite systems written in the left-hand side (31), and the operator fa is determined 
as 



Ta 



■ , 

, 



Substituting Eq. (36) to (34), we have 

^0^7/0 = j d^r Sp2[{tr • r)^fo{6s){r)ra(Ti^il'^(ss){r)] 2i m(£^ - e/J. 
This formula for the dipole emission of photon is similar to that of (13) for the one-channel model. 



(37) 



(38) 



2. Partial width of the decay 4> — »■ 7/0 

Partial width of the decay ^ /o in case, when </> is pure ss state, is determined by the following formula: 

1 — rrvj^ 2 

m^F^^^/o = -a —2 — - \A^^jf^(^ss)\ , (39) 

o rn^ 

where a = 1/137 and the ^^^^^^(55) amplitude is determined by Eq. (6) (here it is specified that we deal with ss 
quarks in the intermediate state). 
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B. Three-component model (ss, nn, gg) for /o and 



The above formula can be easily generalized for the case, when /o is the three-component system (ss, nn, gg) and 
(j) is two-component one (ss, nn), while gg is supposed to be negligibly small. We have two transition form factors: 

F4^^ih{ss) ^ j [(^r)^/„(,,)(r)r„a^V0(.s)(r)] 2z m(e^ - e/J (40) 

and 

p^-.-rMnn) ^ yrf3r-Sp2 [(<Tr)V'/o(„«)(rKa^V0(n^)(r)] 2zm(£^-e/J. (41) 
The partial width reads 

1 m| — nij^ 2 
m^T^^^jg = -a —2 — - + , (42) 

with defined by Eqs. (4) and (6). The charge factors, which were separated in Eq. (4), are equal to: 

y(sg) _ 2 y(nn) _ 1 ..„N 



they include the combinatorics factor 2 related to two diagrams with the photon emmision by quark and antiquark, 
see [16,17] for more detail. 



V. DECAY </-(1020) ^ 7Fo(980) 

The vector meson (/)(1020) has rather small decay width, r0(io2o) — 4.5 MeV; from this point of view there is no 
doubt that treating (/)(1020) as stable particle is reasonable. As to /o(980), the picture is not so determinate. In the 
PDG compilation [20], the /o(980) width is given in the interval 40 < Tfg(^gso) < 100 MeV, and the width uncertainty 
is related not to the data unaccuracy (experimental data are rather good) but a vague definition of the width. 

The mass and width of the resonance are determined by the polo position in the complex- mass plane, M = m—iT/2, 
— just this magnitude is a universal characteristics of the resonance. 



A. The /o(980): position of poles 

The definition of the /o(980) width is aggravated by the KK threshold singularity that leads to the existence of 
two, not one, poles. According to the ii'-matrix analyses [1,5], there are two poles in the {IJ^^ = 00++) wave at 
s ~ 1.0 GeV2, 

~ 1020 - z40 GeV , M" ~ 960 - i200 GeV , (44) 

which arc located on the different complex- Af sheets related to the _ft'_R'-threshold, see Fig. 4. By switching off the 
decay /o(980) KK, both poles begin to move to one another, and they coincide after switching off the KK channel 
completely. Usually, when one discusses the /o(980), the resonance is characterised by the closest pole, M'. However, 
when we are interested in how far from each other the (^(1020) and /o(980) are, one should not forget about the 
second pole. 

Keeping in mind the existence of two poles, one should accept that 0(1020) and /o(980) are considerably "separated" 
from each other, and the /o(980) resonance can hardly be represented as stable particle - we return to this point once 
more in Section 6 discussing the tttt spectrum in (/)(1020) jmr. 
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B. Switching off decay channels: bare states in /f-matrix analysis of the (IJ^'~' = 00"*""*") wave 

A significant trait of the Jr-matrix analysis is that it also gives us, along with the characteristics of real resonances, 
the positions of levels before the onset of the decay channels, i.e. it determines the bare states. In addition, the K- 
matrix analysis allows one to observe the transform of bare states into real resonances. In Fig. 5, one can see such a 
transform of the 00++-amplitude poles by switching off the decays /o — > Tnr,KK,r]r},'r]'r]' ,1^1^. It is seen that, after 
switching off the decay channels, the /o(980) turns into stable state, approximately 300 MeV lower: 

/o(980) f^''"' {700 ± 100). (45) 

The transform of bare states into real resonances can be illustrated by Fig. 6 for the levels in the potential well: 
bare states are the levels in a well with impenetrable wall (Fig. 6a); at the onset of the decay channels (under-barrier 
transitions. Fig. 66) the stable levels transform into real resonance. 

Figure 7 demonstrates the evolution of coupling constants at the onset of the decay channels: following [21], relative 
changes of the coupling constants are shown for /o(980) after switching on/off the decay channels. The onset of the 
decay channels is regulated by the parameter x, and the value x = corresponds to the bare state (amplitude pole 
on the (Re M)-axis) and the value x = 1 stands for the resonance observed experimentally. 

Let us bring the attention to a rapid increase of the coupling constant fo KK on the evolution curve /o^'°(700)- 
/o(980) in the region x ^ 0.8-1.0, where 7^(2; = 1.0) — 7^(a; = 0.8) ~ 0.2, see Fig. 7. Actually this increase allows 
one to estimate a possible admixure of the long-range KK component in the /o(980): it cannot be greater than 20%. 



C. Calculation of the decay amplitude 0(1020) 7/o(980) 

The above discussion of the location of the amplitude poles of /o(980) as well as the movement of poles by switching 
off the decay channels tell us definitely that the smallness of the amplitude of the ^(1020) — *■ 7/o(980) decay due to a 
visible proximity of masses of vector and scalar particles is rather questionable. As to /o(980), its poles "dived" into 
complex plane, in the average in ^100 MeV (40 MeV for one pole and 200 MeV for another). But when we intend 
to represent /o(980) as a stable level, one should bear in mind that the mass of the stable level is below the mass of 
(;6(1020) in 300 MeV — this value is given by the /^-matrix analysis. In both cases we deal with the shifts in mass 
scale of the order of pion mass, that is hardly small in hadronic scale. 

The iiT-matrix amplitude of the 00++ wave reconstructed in [1] gives us the possibility to trace the evolution of the 
transition form factor (I){1020) 7/^'"'^(700 ± 100) during the transformation of the bare state f^^'''{700 ± 100) into 
the /o(980) resonance. Using the diagrammatic language, one can say that the evolution of the form factor P'^^^j^ is 
due to the processes shown in Fig. 8: (j) meson goes into /o"^'^'^(n), with the emission of a photon, then /o'''^°(n) decays 
into mesons fo^^^in) hh = tttt, KK, rjr], rjr]', -k-k-k-k. The decay yields may rescatter thus coming to final states 

The residue of the amplitude pole 0(1020) ^nn gives us the transition amplitude 0(1020) — > 7/o(980). So, in 
the if-matrix representation the amplitude of the reaction 0(1020) 'yrnr, Fig. 8, reads 

p(6are) / s 

A ^ \ 0(lO2O)^7/o''^"(n) bare/ ^ 1 \ /.rn 

Here M„ is the mass of bare state, ga'^'''^{n) is the coupling for the transition f^'^'^^in) a, where a = tttt, KK, r]r], ijr]', 
TTTTTTTr. Thc matrix element {l — ip{s)K{s))~^ takes account of the rescatterings of the formed mesons. Here p{s) is the 
diagonal matrix of phase spaces for hadronic states (for example, for the tttt system it reads pTr-w{s) = {s — 4,m'^)/s), 
and the ii'-matrix elements Kab{s) contain the poles corresponding to the bare states: 

Kab[s) = 2^ -s J'^b[s). (47) 

n " 

The function fab{s) is analytical in the right-hand half-plane of the complex-s plane, at Res > 0, see [1] for more 
detail. 

Near the pole corresponding to /o resonance (resonance poles are contained in the factor (1 — ip{s)K{s))~^ ), the 
amplitude 0(1020) — > 77r7r is written as follows: 

^<?i(i020)^77r7r(s) ^ ^0(10^20)^7/^(980) g ^^^^^^^^^^ smooth terms , (48) 

-'"/o(980) ^ 
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where M/p(980) is the complex-valued resonance mass: M^j^fggo) ~ 1020 — i40 MeV for the first pole, and 

■^/o(980) ~^ — 960 — «200 MeV for the second one. The transition amplitude ^0(1020)^7/0(980) is different for 
different poles; the 5/(,(980)^7r7r couplings are different as well. 

We see that the radiative transition (^(1020) 7/0 (980) is determined by two amplitudes, 
^0(iO2O)^7/o(Mi) = ^0(1020)^7/0(980) ^nd ^0(io2o)^7/o(M") = ^0(io2o)^7/o(98O) ' and just thesc amplitudes are the 
subject of our interest. The amplitudes ^0(1020)^7/0(980)' "^0(1020)^7/0(980) '^ay be represented as the sum of contri- 
butions from different bare states: 

^0(1020)^7/0(980) = X]^"[-^o(9^^)]-^0(lO2O)^7/ob"-(n)' (^9) 
n 



^<?!>(1020)^7/o(980) ~ 



To calculate the constants CnUoimR)] we use the ii'-matrix solution for the 00++ wave amplitude denoted in [1] as 
II-2. In this solution, there are five bare states fg^'^'^in) in the mass interval 290-1950 MeV: four of them are members 
of the qq nonets {l^Poqq and 2^Poqq) and the fifth state is the glueball. Namely: 



2''Paqq 
glueball 



/o'''"'°(700± 100), 
/□'"^'■°(1230±40), 
/(5"^'°(1580±50). 



/(5"'''*'(1220±30), 
/o'"'''''(1800±40), 



(50) 



For the first pole of the /o(980) resonance located at M[/o(980)] 
as follows: 



1020 — i40 MeV the renormalization constants are 



41,1, [/o (980)] 
Cf2^2o[/o(980)] 
d2'3o[/o(980)] 



0.62exp(-a44°), 

0.37exp(-i41°), 

0.19exp(a°), 



(51) 



Co[/o(980)] = 0.02exp(-ai2°), 
Cf58o[/o(980)] = 0.02exp(i5°). 

These constants are complex- valued. One should pay attention to the fact that the phases of constants C7oo[/o(980)] 
and Ci22o[/o(980)] have the relative shift close to 90°. This means that the contributions from /q'''''^(700 ± 100) and 
jbarc 1^2220 ± 30) (which are members of the basic l^Poqq nonet) do not interfere practically in the calculation of 
probability for the decay ^(1020) ^ 7/o(980). 

Actually one may neglect the bare states f^'^'^{1230), /^"""^{ISOO), f^'^'^{1580) in the calculation of the (?i(1020) ^ 
7/0 (980) reaction, bacause the form factors for the production of radial excited states are noticeably suppressed, see 
[17]: 



p(6are) 

■ 0(lO2O)^7/o(23Pogg) 



F 



(bare) 



0(1020)^7/0 (l^'i'oqq) 



Besides, the coefficients d23o[/o(980)] d8oo[/o(980)] are also comparatively small, see (??). 

The second pole located on the third sheet, M[/o(980)] = 960 — z200MeV, has renormalizing constants as follows: 

4|,'J[/o(980)] = 1.00exp(26°), 
d22o[/o(980)]=0.33exp(zll3°), 
d23o[/o(980)]=0.32exp(a48°), 
d80o[/o(980)] = 0.08exp(i4°), 
Cf58o[/o(980)] = 0.04exp(i98°). 

Here, as before, the transitions (/)(1020) ^ 7/o''"°(1230), 7/o'"'''''(1580), 7/o^^''°(1800) are negligibly small. 



(52) 
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In (/)(1020), the admixture of the nn-component is small. In the estimates given below we assume 0(1020) to be 
pure ss state. The bare states fQ^'"^{700) and /q '^''®(1220) are mixtures of the nn and ss components 

nn cos (fi + sssinip, 

and, according to [1], the mixing angles are as follows: 

^ [/bare (700)] = -70° ± 10°, 

if [f^'''%1220)] =20° ±10°. (53) 
Because of that the transition amplitude for 0(1020) — > 7/o(980) reads 

^^(1020)^7/0(980) - Ctoo[./o(980)] sin [/(?"°(700)] ^^(i020)^7/b»-(700) 

+ Cf22o[/o(980)] sin^ [/^^'■^(1220)] i^^;:;2i)^7/--(i220) • (54) 

Here C7oo[/o(980)] and Ci22o[/o(98G)] {N = I, II) are given by the formulae (??), (52). One can see that numerically 
the factor Ci22o[/o(980)] sini^ [/o''™(1220)] is small, and we may neglect the second term in the right-hand side (54). 
Then for the pole, which is the closest one to the real axis (1020 - 140 MeV), one has: 

^Ji(1020)^7/o(980) - (0-58 ± 0-04)i^^("o20)^7/b="-''(700) ' (^5) 

and for the distant one, (960 - 1200 MeV): 

4(1020)^7/0(980) - (0.92 ± 0.06)f(;^;^)^^^^,.„(,„„) . (56) 

We see that practically the ^^^(1020)^7/0(980) amplitude does not change its value during the evolution from bare state 
to resonance, while the decrease of ^^1(1020)^7/0 (980) significant. 

D. Comparison to data 

Comparing the above-written formulae to experimental data we have parametrized the wave functions of the qq 
states in the simplest, exponent-type, form (see Section 2.3). For 0(1020), we accept its mean radius square to be 
close to the pion radius, ^^(1020) — -^w- both states are members of the same 36-plet. This value of the mean radius 
square for 0(1020) fixes its wave function by = 10 GeV~^. 
For /o^''''(700), we change the value in the interval 



5GeV-'<6(/"''^)<15GeV- 

— JO — 



that corresponds to the interval (0.5-1.5)i?J of the mean radius square of Z^""''^) (700). 
Using the branching ratios [22,23] as follows: 



BR[0(1O2O) ^ 7/0 (980)] = (3.5 ± 0.3t J j) x 10"^ , (57) 

BR[0(1O2O) ^ 7/o(980)] = (2.90 ± 0.21=^^-^^) x 10"^ , 
and the definition of the radiative decay width: 

1 m| — mj^ 2 
m^r^^^/o = —2 1^0^7/0 1 , 

we have the following experimental value for the decay amplitude: 

^£7o20)^7/o(980) = 0.115 ± 0.040 GeV . (58) 
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Here a = 1/137, = 1.02 GeV and mj„ = 0.975 GcV (the mass reported in [22,23] for the measureed 7/o(980) 
signal) and rtot[0(lO2O)] = 4.26±0.05 MeV [20]. The right-hand side (58) should be compared with ^^(1020)^7/0(980) 
calculated with Eqs. (17), (38), and (55): 

4ri020)^-y/o(980)(dipole) (0.58 ± 0.04) ^W^,,[/b-'^ (700)] Z^^^ x (59) 



^7r(6;t^j^-^[-^-(«-^^°-^)^^^]- 

Recall, in (59) the factor (0.58 ± 0.04) takes into account the change of the transition amplitude caused by the final- 
state hadron interaction, Eq. (55). The probability to find quark-antiquark component in the bare state /o'''^''(700) 
is denoted as Wqq[fQ^'"^ (700)]: one can guess that it is of the order of 80-90%, or even more. The mass of the strange 
constituent quark is equal to nis — 0.5 GeV. The wave functions of ^(1020) and /o'^'''^(700) are parametrized as 
exponents: we fix = 10 GeV~^ (that gives for the mean radius of 0(1020) the value of the order of the pion radius 
, i?0 ~ Rt:), and vary bf„ in the interval (5-15) GcV~^. 

The comparison of the data (58) to the calculated amplitude is shown in Fig. 9. We see that the calculated 

amplitude (59) is in a perfect agreement with data, when M,ibare) ~750-800 MeV, that is just inside the error bars 

•'0 

given by the JC-matrix analysis [1]. 



VI. PION-PION SPECTRUM IN 0(1020) -> itttt 



The /o(980) resonance is seen in the reaction 0(1020) — > 77r7r as a peak at the edge of the tttt spectrum. So it 

is rather enlightening to calculate the tttt spectrum to be sure that its description agrees both with the quark model 



■o(980) s-nd threshold theorem (cross section tending to zero as w at a; ^ 0, 



calculation of the form factor ^0(1020)^7/1 
where lo = — M^jr)- 

Partial cross section of the decay 0(1020) 777^77° is given by the following formula: 



dr. 



<i>{W20)^1vOw° _ 1-p '•^'^ - /l/r2 
- 0(1020)^7/0(980) —2 



"''<P ^'-^TTTT 



(60) 



Jo 



2M^ 



97r 



The factor 1/3 in front of the right-hand side (60) is associated with the tt^tt^ channel: r0(io2o)^77rO7r'' = 
l/3r0(io2o)^77i-7r- Here for the description of the /o(980) we use the Flatte formula [24] with the phase space factors 



(61) 



At M^^ < 4m|^ one should replace M^^ - 4to|. i^Am\ - M^^. In line with [22,23,25], we use the Flatte 
formula with the parameters 



gl = 0.12 GeV^ , g% = 0.27 GeV^ , Mq = 0.975 GeV . 

The threshold theorem requires 

.97r 



M2 - Ml 



+ B{M, 



tttt) 



(M^„ - m^) 



(62) 



(63) 



that gives a constraint for the background term Bi^M"^^). The term Bi^M"^^) is parametrized in the form: 



B{Ml,)^C[\ + a{M., 



2 

TTTT 



ml)] exp 



(64) 
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and the parameter C is fixed by the constraint: 



djL + (7 



= . (65) 



Fitting to the -k^-k^ spectrum [22], see Fig. 10a, we have the following values for other parameters: 

^ = -0.2GcV^ = 0.388 GcV . (66) 

For r0(io2O)^7/o(98O) entering (60) we have used ^,^(io20)-»7/o(980) = 0-13 GeV that satisfies both (58) and (59). 

The Flatte formula gives us rather rough description of the tttt amplitude around the /o(980) resonance. More 
precise description may be obtained by using in addition the non-zero transition length for tttt — > KK [21]. For this 
case, we have the formulae analogous to Eq. (60), after replacing the resonance factor 

M2 - - iglp^^ ~ igjiPKK ^^^^ 

by the following one: 

+ iPKRdKf 



The parameters found in [21] are equal to: 



(68) 



= 0.386 GeV, gx = 0.447 GeV, (69) 
Mo = 0.975 GeV,/ = 0.516. 

The transition length a^^^j^j^ is determined by the parameter / as follows: (It^t^^kk = '^f l^o- 

The description of the 7r°7r° spectra [22] within the resonance formulae (68) is demonstrated in Fig. 106. In this fit 
we have the following parameters for B{M'^^): 

a = 0, A* = 0.507 GeV. (70) 

In this variant of the fitting to spectra wc also used ^0(1020)^7/0(980) = 0.13 GcV . 

Let us emphasize that the visible width of the /o(980) signal in the tttt spectrum is comparatively large, ~150 MeV, 
that is related to an essential contribution of the second pole at 960 — i200 MeV. 



VII. THE ADDITIVE QUARK MODEL, DOES IT WORK? 



Let us point to the two aspects of this question. One is the problem of the applicability of the additive quark 
model to the production of the resonance /o(980), another one is the production of the bare state /q (700 ± 100). 



A. Process 0(1020) ^ 7/i""'''(700 ± 100) 

The additive quark model describes well the production of the bare state /o**"^® (700 ±100), provided its mass is in the 
region 750-800 MeV. To see it, consider Eq. (38) for F^^''^" (dipole), or Eq. (17), where exponential representation 
of the quark wave functions is used. Formula (17) takes into account both the additive quark model processes and 
photon emission by the charge-exchange current, while Eq. (16) gives us the triangle-diagram contribution within 
additive quark model. The contribution of the charge-exchange current is small, when 

ms[m^ - Afybare] ~ (71) 

At TTig = 0.5 GcV and = 10 GcV~^ the equality (71) is almost fulfilled, when M^bare ~ 0.8 GeV. Such a magnitude 
is allowed by the if -matrix fit [1], which gives Mjbare = 0.7 ± 0.1 GeV. 

However, let us emphasize that the error bars ±0.1 GeV are rather large in the difference (m^ — M^bare): with the 
lower possible limit Mjbare = 0.6 GeV we face a two-times disagreement in Eq. (71). Still, one may hardly hope that 
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the ii'-matrix analysis of the 00"*"^ wave would provide us with a tighter restriction for the mass of this bare state, 
since a large uncertainty in the definition of Mjbare is not related to the data accuracy but to the problem of the light 
a meson existence, see the discussions in [5,26-28] and references therein. 

The use of i^^^ (additive), Eq. (16), for the calculation of ^^^(^ 020) ^7/0(980) I'^sults in the agreement with 
experimental data. Thus we have: 

<(i020)^,/„(980) (additive) ^ (0.58 ± 0.04)^^^^^^] 4^7/o (^2) 



2V^ 6f5f 



%/3 (6^ + 6/o)^/^' 

In Fig. 11, one can see ^^11020)^7/0(980) (additive) versus ^^'(^020)^7/0(980) ' there is a good agreement with data. 

The existence of two characteristic sizes in a hadron, namely, hadronic radius and that of constituent quark, may 
be the reason why the contribution of the charge-exchange current is small in the reaction 0(1020) 7/0''''® (700). 
Relatively small radius of the constituent quark assumes that charge-exchange interaction ss ^ gg nn is a short- 
range one, that causes a smallness of the second term in the right-hand side (35). 

The hadronic size is defined by the confinement radius ~ -Rconf , which is of the order of 1 fm for light hadrons. 
The constituent quark size, r^, is much smaller, it is defined, as one may believe, by relatively large mass of the soft 
gluon (experimental data [29] and lattice calculations [30] give us nig ^700-1000 MeV). So we got r^/R^ ^(0.1-0.2), 
the same value follows from the analysis of soft hadron collisions, see [15,31] and references therein. 



B. Process </)(1020) 7/o(980) 

The two sizes, and i?^, being accepted, the additive quark model contribution dominates the reaction (^(1020) 
7/0 (980) too, thus allowing direct use of the triangle diagram of Fig. la for the calculation of this process. Such 
calculations were performed in [17], revealing resonable agrement with data. Once again it should be emphasized that 
the triange diagram contribution does not have a particular smallness related to a deceptive proximity of (j){1020) and 
/o(980). Besides, as was explained above, the poles associated with these resonances are separated from each other 
in the complex- Af plane in nonsmall distances in the hadronic scale. 

In the literature there exist rather opposite statements about the possibility to describe the reaction 0(1020) 
7/0 (980) within the frame of the hypothesis of the qq nature of /o(980). Using the QCD sum-rule technique the 
authors of [32] evaluated the rate of the decay 0(1020) — > 7/o(980), with a fair agreement with data, supposing a 
sizeable ss component in the /o(980). 

The results of the calculation performed in [33] in the framework of the additive quark model do not agree with 
data on the reaction 0(1020) 7/0 (980). This calculation though similar to those of [16,17] led to different result, 
so it would be instructive to compare model parameters used in these two approaches. 

In [33] as well as [16,17], the exponential parametrization of the wave function was used, however the slopes 
and bfg in [33] were considerably smaller (constituent quark masses are smaller too). In [33], bss — 2.9 GeV~^ and 
= b^d 3.7 GoV-2 (m„ = nid = 220 MeV, m., = 450 MeV), while in [16,17] 6^ ~ 10 GeV'^ and b^, - b^ 
(m„ = rrid = 350 MeV, rUg = 500 MeV). Besides, in [33] the scheme of the mixing of /o states was used that was 
suggested in [34,35], where the transitions Jq^'^'^ real mesons were not accounted for. Still, as was emphasized above 
(Section 5.2), just the transitions /q""'® — > tttt, KK, -qr], TrTTTTTr afford the final disposition of poles in the complex plane, 
for they are responsible for the resonance mass shift of the order of 100 MeV, see Fig. 5. 

In our opinion, the failure of the qq model demonstrated in [33] can testify only the fact that not any model enables 
the description of radiative decays. The qq model should be based on the whole set of experimental data but not on 
the reproducing several levels of the lowest states. 



VIII. CONCLUSION 



Correct determination of the origin of /o(980) is a key for understanding of the status of the light a and classification 
of heavier mesons /o(1300), /o(1500), /o(1750) and the broad state /o(1200-1600). 

We have shown that experimental data on the reaction 0(1020) 7/o(980) do not contradict the suggestion about 
the dominance of the quark-antiquark component in the /o(980). However, as was emphasized in Introduction, the 



15 



final conclusion about the origin of /o(980) should be made on the basis of the whole availability of arguments, so let 

us enumerate them briefly. 

1. There are data on the hadronic decays of the lowest mesons, and the most reliable information on scalar 
resonances is given by the if -matrix analysis. Summing up, one can state that, according to the if -matrix analysis 
[1], the lowest states /o(980), /o(1300), ao(980), isro(1430) are the descendants of bare mesons, which have created 
the l^Po multiplet. Because of that, all the decays, namely, 

/o(980)^7r7r,if^, (73) 
/o(1300)^7r7r,if^,W, 

ao(980) KK,n'n, 
ifo(1430) isTTT 

are described by two parameters only in the leading terms of the 1/Nc expansion, which are the universal coupling 
constant and mixing angle for the nn and ,ss-components in the scalar-isoscalar sector. 

2. Another argument is the systematics of scalar states on linear trajectories on the (n, Af^) plane. For scalar 
mesons the trajectories are shown in Fig. 12a — one can see that /o(980) lays comfortably on linear trajectory, 
together with the other scalars. Such trajectories are formed not only for the scalar sector but also for all aggregate 
of data, see [4,5], and all the trajec;tories are characterized by the universal slope. 

Similar qq trajectories exist in the (J, M^)-plane too, and /o(980) belongs to one of them. 

3. The alternative to the qq system may be the four-quark qqqq, or molecular KK structure [6,7,36,37]. Such a 
nature of /o(980) would mean that /o(980) was loosely bound system but the experiment tells us that this is not so. 
The matter is that 

(i) /o(980) is easily produced at large momenta transferred to the nucleon in the reaction 7r~p — > /o(980)n [38,39], 

(ii) /o(980) is produced in the Z^-boson decays [40], 

(iii) /o (980) is produced in central pp collisions at high energies [41] . 

Were the /o(980) a loosely bound system, these processes would be suppressed. 

4. The /o(980) resonance is produced in hadronic decays of the meson, i?+ 7r+/o(980), with the probability 
comparable with that for the transition £)+ 7r"'"(?!i(1020) [42]. These two reactions are due to the weak decay of 
c quark, c tt+.s: as to /o(980), it is formed, like (^(1020), by the ss system in the transition ss /o(980). The 
calculation of this process [43] shows us that the /o(980) yield in the reaction i)+ 7r+/o(980) can be reliably 
calculated under the assumption that /o(980) is close to the qq flavor octet. 

The study of the Df 7r+/o(980) decay by [44-46] also led to the conclusion about the ss nature of /o(980). 

5. Concerning radiative decays with the formation of /o(980), we see that the transition 0(1020) 7/o(980) can 
be well described within the approach of additive quark model, with the dominant qq component in the /o(980). 
Another radiative decay, /o(980) 77, partial width of which was measured [47], can be also treated in terms of the 
qq structure of /o(980) [16,48]. The values of partial widths in both decays support the conclusion made in [1] that 
the flavor content of /o(980) is close the octet one. 

We thank Ya.I. Azimov, L.G. Dakhno, S.S. Gershtein, Yu.S. Kalashnikova, A.K. Likhoded, M.A. Matveev, V.A. 
Markov, D.I. Melikhov, A.V. Sarantsev and W.B. von Schlippe for helpful and stimulating discussions of problems 
involved. 

This work is supported by the RFBR Grant N 04-02-17091. 



APPENDIX A: DIPOLE EMISSION OF PHOTON 



To describe the interaction of composite system with electromagnectic fleld, we consider the full Hamiltonian: 



ii(0) 



2m ' 2m 



Uss^ssiri - r2) , 



Us. 



?(ri - ^2) 



ss—ygg 



(ri-r2) , it + ^ + U.. 



(74) 



Here the coordinates (r^) and momenta (ka = — zVa) of the constituents are related to the characteristics of the 
relative movement, entering (31), as follows 



ri = ir-t-R, r2 = -ir-FR, ki = ik-hP, k2 = -ik-FP. 



(75) 



The electromagnetic interaction is included by substituting in (74) as follows: 
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(ki - eiA(ri)) 



(k2 - e2A(r2))' 



(76) 



- r^) Uss^gg{vi - T2) exp 



ri r2 

iei j dr'^Aair') + ie2 j dr'^A^ir') 



with ei = — 62 = eg. After that we obtain the gauge- invariant Hamiltonian H{A): 

H{A) = x+H{A + Vx)x, 
where A + Vx means the following substitution: 

A{ra) A{ra) + Vx(ra) , 

and matrix x reads 



X 



, 1 



(77) 
(78) 

(79) 



For the transition — > 7/0, keeping the terms proportinal to the s-quark charge, eg, we have the following operator 
for the dipole emission: 



2(fcia - k2a) , i{ria - ria)Uss^gg{ri - r2) 

-i{ria - ria)Uss^gg{ri - ra) , 



(80) 



There exist other mechanisms of the photon emission which, being beyond the additive qiiark model, lead lis to the 
dipole formula for V transition, an example is given by (L • S)-interaction in the quark antiquark component 

[9,10,49]. The short-range (L • S) interaction in the qq systems, was discussed in [50,51] as a source of the nonet 
splitting. Actually the point-like (L • S) interaction gives {v/c) corrections to the nonrelativistic approach. In the 
relativistic quark model approaches based on the Bethe-Salpeter equation, the gluon-exchange forces result in similar 
nonet splitting as for the (L • S) interaction, for example, see [52] . 
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FIG. 3. Examples of diagrams, which contribute to the potentials Uss^ss{r), Uss^gg{r) and Ugg^gg{r). 
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FIG. 4. Complex-M plane and location of the poles corresponding to /o(980); the cut related to the KK threshold is shown 
as broken line. 
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FIG. 5. Complex M plane: trajectories of poles corresponding to the states /o(980), /o(1300), /o(1500), /o(1750), 
/o(1200 — 1600) within a uniform onset of the decay channels. 
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FIG. 6. The /o-levels in the potential well depending on the onset of the decay channels: bare states (a) and real resonances 



(b). 
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FIG. 7. The evolution of normalized coupling constants 7a = ffa/y^;, gl at the onset of the decay channels for /o(980). 
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FIG. 8. Diagram for the transition ^(1020) — > ^tttt in the ii'-matrix representation, Eq. (46). 
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FIG. 10. The tttt spectrum of the reaction 0(1020) '^tx'k calculated with the Flattee formula (a) and Eq. (68) (b). 
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FIG. 12. Lineax trajectories on the (n, M^) plane for bare states (a) and scalar resonances (b). 
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